Orbital Selectivity in Hund's metals: The Iron Chalcogenides 
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We show that electron correlations lead to a bad metallic state in chalcogenides FeSe and FeTe 
despite the intermediate value of the Hubbard repulsion U and Hund's rule coupling J. The evolution 
of the quasi particle weight Z as a function of the interaction terms reveals a clear crossover at U ~ 
2.5 eV. In the weak coupling limit Z decreases for all correlated d orbitals as a function of U 
and beyond the crossover coupling they become weakly dependent on U while strongly depend on 
J. A marked orbital dependence of the Z's emerges even if in general the orbital-selective Mott 
transition only occurs for relatively large values of U. This two-stage reduction of the quasi particle 
coherence due to the combined effect of Hubbard U and the Hund's J, suggests that the iron-based 
superconductors can be referred to as Hund's correlated metals. 

PACS numbers: 71.30.+h, 71.10.Fd, 71.27. +a 



The role of electron correlations in the iron-based su- 
perconductors is still a debated issue, naturally inter- 
twined with the search for the origin of high critical tem- 
peratures. We present results that improve the qualita- 
tive understanding of how electron correlation influences 
fundamental electron properties of these compounds, 
such as the metallicity, which in turn might be impor- 
tant also for the understanding of the pairing mechanism. 
We choose two candidates of the chalogenides, FeSe and 
FeTe and employ first principles electron structure cal- 
culations combined with advanced many-body methods 
taking into account the local electron correlation. The 
chalcognides have in contrast to the pnictides a simpler 
atomic structure, thus easier to synthesize and also to 
study theoretically. In addition they are non toxic in 
contrast to the pnictides containing arsenic. 

In previously known superconductors we can identify 
either weakly correlated materials, like elemental super- 
conductors or binary alloys, including MgB2, or highly- 
correlated compound like the copper oxides and heavy 
fcrmion materials. In the first set of compounds su- 
perconductivity is explained within the Bardeen-Cooper- 
Schrieffcr framework and its extensions, and it occurs as 
a pairing instability of a normal metal. In the second set 
it is widely believed that correlations revolutionize the 
electronic properties and that both the metallic state and 
the pairing mechanism deviate from standard paradigms. 

The iron-based pnictides and chalcognides supercon- 
ductors do not fit this simple classification. The common 
labeling "intermediate correlation" , referring to proper- 
ties such as Fermi surface topology or absence of Hub- 
bard bands [l|, suggests modest effects of correlations. 
Conversely, the metallic state appears much less coherent 
than what these observations imply 0, E| . A magnetic 
counterpart of this dualism is the localized an itinerant 



nature of the spin-density- wave state of the parent com- 
pound. 

The characteristic property of the band structure is 
that several of the five d-bands cross the Fermi level. 
The multi-orbital nature leads to several exotic electronic 
properties such as orbital-selectivity 043 an d a ls° to the 
conclusion that the inter-orbital exchange or Hund's cou- 
pling plays a key role 13, 

The role of the Hund's coupling has indeed been recog- 
nized in the early stages of the field in a Dynamical Mean- 
Field Theory (DMFT) study by Haule and Kotliar [nj, 
who coined the definition of Hund's metals by the ob- 
servation that the quasi particle effective mass and the 
response functions are much more sensitive to the Hund's 
coupling J than to the Hubbard U interaction. 

For a Hund's metal the spectral weight is not trans- 
ferred to the high-energy Hubbard bands, but rather 
spreads over a scale controlled by J. Other DMFT stud- 
ies have highlighted the anomalies of the metallic state, 
showing its incoherent nature [l3j and its relation 
with a spin-freezing crossover jl4 |. In Ref. [l5[ the dual 
nature of the magnetic correlation is shown to induce 
a remarkable difference between a large instantaneous 
magnetic moment and smaller long-time magnetic cor- 
relations, similar to the spin-freezing scenario proposed 
in Ref. [l6| for a three-orbital model. 

It has recently been shown that J can have a two-fold 
effect on a multiorbital system with an integer filling dif- 
ferent from one electron per orbital [I?} , a situation which 
is realized in the parent compounds of iron superconduc- 
tors, in which six electrons populate the five d orbitals. 
In this configuration J reduces the quasi particle coher- 
ence temperature (or coherence energy scale), while it 
increases the critical U for the Mott transition. 

As a consequence, a two-stage reduction of the elec- 
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tronic coherence scale (measured by the quasi particle 
weight Z) occurs as a function of U. Indeed, if we choose 
a sizable value of J and follow the evolution of the metal- 
lic properties, we first have a rapid decrease of the effec- 
tive Fermi-liquid coherence scale, which leads to a bad 
metal already for intermediate correlations strengths, 
while the Mott transition occurs only at much larger U. 
This opens a window of U in which Z is essentially flat, 
which has been dubbed after the roman god Janus in view 
of the double-faced effect of the Hund's coupling 17 1. 

In this work we explore the combined role of U and J 
in the iron-based chalcogcnidcs FeSe and FeTe by means 
of the Gutz wilier approximation (GA). The GA is a sim- 
plified treatment of electron correlations which system- 
atically selects the energetically favorable electronic con- 
figurations out of an uncorrelated wave function. The 
method provides a reasonable description of the Mott 
transition from the metallic side [l8[ and allows for a nu- 
merically cheap investigation of a wide range of model 
parameters. 

We employ the GA numerical scheme developed in 
Ref. [19l42l| , which is a generalization of earlier formula- 
tions of the GA method [22M25I] and which enables tak- 
ing into account the full rotationally invariant Hund's 
terms, including the so-called spin-flip and pair-hopping, 
that are often hard to treat with approximate analytical 
methods and even with numerical methods. Since the 
formation of a Hund's metal is actually associated with 
a differentiation between the different atomic multiplets, 
we expect that the GA will perform even better than for 
standard Mott transitions. 

Based on electronic structure calculations of FeSe and 
FeTe combined with the GA we show that the electronic 
configuration of the parent compounds of iron-based su- 
perconductors form an ideal system with a two-stage re- 
duction of electronic coherence. Furthermore, the bad 
metal arising from the interplay of U and J displays, as 
expected, an orbital-selective coherence with t2 g orbitals 
significantly more correlated than e g . 

The material-specific band structure is determined us- 
ing Density Functional Theory with the Generalized Gra- 
dient Approximation for the exchange-correlation poten- 
tial according to the Perdew-Burke-Ernzerhof recipe as 
implemented in Quantum Espresso [26| . Then we apply 
Wannier90 |27J to compute the maximally localized Wan- 
nier orbitals, and we include the interaction terms of the 
form 

H = U ^2 nimfriimi + (V - Tj-) n im n im , (1) 

4,+ H.c.) . 

i,m>m' 

Here di^ ma is the destruction operator of an electron of 
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is the spin operator for orbital 
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FIG. 1: Quasiparticle weight in a N or b-fold degenerate Hub- 
bard model with average population of N or b+l electrons, as 
a function of U/D for various Hund's coupling J/U, where 
D is the half-bandwidth. Upper panel: N or b=5, Lower panel 
N or b=3. The panels show data for N=N or b+l electrons, the 
insets for N=N orb . 



m at site i. U and U' = U—2J are intra- and inter-orbital 
repulsions and J is the Hund's coupling. The values of U 
and J are not directly accessible from experiments and 
even if reliable theoretical estimates are obtained with 
constraincd-RPA, there are still some discrepancies be- 
tween different calculations. In the light of the extreme 
sensitivity on the value of the parameter J, it is partic- 
ularly useful to apply a method such as the present GA 
which allows for a continuous sweep of many parameter 
values. 

As mentioned above, one of the distinctive features of 
the iron-based superconductors is that all five d-orbitals 
appear to contribute to the band structure around the 
Fermi level. As a first step we consider a system with 
five degenerate d bands, and show that the configuration 
with six electrons per atom, characteristic of the parent 
compounds, is a clear cut case of a "Janus" scenario, 
characterized by a two-stage reduction of the quasi par- 
ticle weight. 

In Refs. [ll, 17 1 it is clearly shown that the two-stage 
reduction of the electronic coherence scale is a conse- 
quence of a contrasting effect of J on the metallic char- 
acter of the electrons. In the weak-coupling limit J favors 
the formation of a large local magnetic moment, which 
leads to a faster decay of the electronic coherence scale 
Z, while in the strong-coupling the Mott transition is 
pushed to larger U . This effect is particularly strong 
when the number of electrons per atom differs by one 
unit from the number of orbitals N = A or b ± 1, and it is 
expected to be emphasized increasing the number of or- 
bitals as the weak-coupling coherence temperature scales 
exponentially with N m \y. 
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FIG. 2: Upper panel: quasiparticle weight in a Hubbard 
model with 6 electrons in 5 bands with semi-circular den- 
sities with half-bandwidth D split by a cubic crystal-field in 
two manifolds of degeneracy 3 (t2 g symmetry) and 2 (e g sym- 
metry). Lower panel: populations of the two manifolds. 



In Fig. Q] we compare the GA results for the two cases 
of N = N OI h ± 1 when iV or b = 5 and 7V or b = 3. We clearly 
see that the former case has a much clearer separation be- 
tween a regime in which Z rapidly decreases as a function 
of U and a large bad metal region in which Z is essentially 
constant prior to the Mott insulator transition. In the in- 
set we show the case of half- filling, N = iV or b, where no 
dual nature is observed. 

Once established that the electron count of the par- 
ent compounds of the iron-based superconductors gives 
rise to a strongly two-faced correlation physics, we move 
towards the realistic situation in order to identify how 
the material-specific properties influence this picture. As 
an intermediate step, we lift the degeneracy with a cubic 
crystal-field which is present in iron pnictides and chalco- 
genides. An energy splitting A is introduced between the 
three ti g and the two e g orbitals. In Fig. [2] we show the 
results for A/D = 0.2. We observe that, while the weak- 
coupling region gives an essentially orbital-independent 
Z, as soon as we enter in the strongly correlated re- 
gion, the low-lying states become more correlated than 
the higher-lying. In other words, the crystal-field trig- 
gers a strongly orbital-selective renormalization in the 
bad metal state. 

We finally perform the realistic DFT+GA calculation 
for iron chalcogcnides. In panel (a) of Fig. [3] we show the 
evolution of the quasi particle weight for the different or- 
bitals function of U, keeping the ratio J/U fixed 
to 0.224. This ratio is chosen according to the estimates 
presented in Ref. [28| for FeSe. The picture remains simi- 
lar to the idealized systems. For small values of U the Z's 
for the different orbitals are similar and they appreciably 
decrease as a function of U before U ~ 2.5 eV, where 
the system enters the novel regime in which the Z's are 
small and almost constant as a function of U. However, 
an orbital dependence also appears clearly. In addition 
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FIG. 3: Results for FeSe (DFT+GA with J/U = 0.224) as 
a function of U. From the top; panel (a), the quasi particle 
weights for the different orbitals, in panel (b) the expectation 
value of S 2 , in panel (c) the renormalized crystal-field split- 
tings, in panel (d) the population of each orbital, in panel (e) 
the inter-orbital density correlations. 



to the differentiation of the ti g and e g orbitals, we find 
that the d xy orbital is the most correlated and the d x 2_ y 2 



is more localized than the d 3z 2_ r 2. The crossover, which 
roughly separates a weakly-correlated phase from a bad 
metallic phase, takes place at a value of U smaller than 
the bandwidth 2D (~ 4eV), and much smaller than the 
multiband Mott transition, which would take place at a 
U of the order of 5 times the width of each band [29| . 

It is easy to see that in the atomic limit the ground 
state changes from low-spin (S = 1) to high-spin (S = 2) 
when J becomes larger than the crystal-field splittings 
(~ 0.6eV for FeSe). In the metallic phase this evolution 
yields a crossover to a S = 2 state (Panel (b)) which 
leads to the rapid reduction of Z [llj . 

This crossover leads to a dramatic lowering of the co- 
herence temperature, and opens a wide bad-metal re- 
gion due to the increase of U c induced by the effect 
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FIG. 4: Results for FeSe and smaller J/U = 0.15 (top panel) 
and for FeTe (bottom panel) 



of J on the high-spin Mott gap. This behavior is ob- 
served even more pronounced in studies of LaFeAsO 

It contributes in 



and the intercalated chalcogcnides [31 1 
a substantial manner to the sharp onset of the Hund's 
metal (Io| / spin- frozen [3] /incoherent [H, 13| phase ob- 
served in all DMFT studies. 

This effect is also reflected in the rcnormalized orbital 
energies (panel (c) of Fig. [3]), with four of the five or- 
bitals being brought close to one another and also near 
the Femi level by the interactions. This favors a more 
even population of the orbitals that gains in exchange 
energy, and favors the high-spin configurations over the 
low-spin ones. The evolution of the population of the 
different orbitals is shown in panel (d) of Fig. [3] The 
t-2 Q orbitals have populations closer to half-filling already 
at the DFT level, while the e g bands are more occupied. 
Increasing the interaction the d xy level, which has the 
smaller Z , becomes less occupied than d xz and d yz due 
to the stronger effect of correlations. Analogously, the 
large difference between the non-interacting densities of 
the two e g bands is washed out by correlations, which 
favor a more democratic occupation with high spin. The 
same low to high spin transition within the metallic phase 
occurs for any sizable value of J and it is indeed present 
already in the model with a simple ti g -e g splitting as 
clear from the non-monotonic population behavior shown 
in the lower panel of Fig. [2j 

In Panel (e) of Fig. [3] we show the inter orbital den- 
sity correlation functions, which are clearly suppressed 
in the correlated regime. This suppression, driven by J, 
has been put forth 0, [2^ as the driving mechanism be- 
hind the orbital selectivity. Indeed J acts as an "orbital 
decoupler" ("band decoupler" (29|, for weak orbital hy- 
bridization) suppressing inter-orbital charge fluctuations, 
and rendering the charge dynamics of each orbital virtu- 
ally independent. For a smaller value of J/U, 0.15, the 
picture does not change. The position of the crossover is 
only weakly affected, while the values of the Z'a in the 



bad metallic region after the crossover depend strongly 
on J. The behavior observed confirms previous findings 
of a Z which depends strongly on J and weakly on U in 
the physically relevant region of U ~ 4 eV and J ~ 0.5 — 1 
cV. At the same time our results clearly underline that 
such "Hund's metal" requires a critical value of the Hub- 
bard repulsion, albeit much smaller than one might ex- 
pect on the basis of the value of the bandwidth. The 
picture is clearly consistent with that drawn in Ref . [13] ■ 
Finally, we show the quasi particle weights for FeTe using 
the same value of the interaction coefficients. The main 
difference is a sharper separation between t2 g and e g or- 
bitals, and a larger renormalization for the e g orbitals. 

In summary, we have calculated the correlation 
strength induced by many-body correlations on the ab- 
initio electronic structure of the iron chalcogcnides FeSe 
and FeTe. We find, in agreement with previous analo- 
gous studies on LaFeAsO and Ki_ a: Fe2-ySe2 [3l| . 
that Hund's coupling has a strong influence on the elec- 
tronic properties of the paramagnetic phase, inducing a 
two-stage quasi particle renormalization. A first regime 
at weak coupling sees a moderate correlation affecting 
all orbitals comparably. After a quick decrease around 
U ~ 2.5eV, a value much smaller than the overall band- 
width, a strongly correlated regime is entered, heavily 
differentiated among the orbitals (with t2 g orbitals sensi- 
bly more correlated), in which the quasi particle weights 
are almost independent of U. The Mott transition occurs 
at much higher (~ 5 times the bandwidth) interaction 
strengths. Comparison with idealized models shows that 
the two-staged reduction of the quasi particle weights is 
due to the filling of 6 electrons in 5 bands, thus placing 
the system in the "Janus" regime induced by Hund's cou- 
pling [17|, and that, by introducing a crystal-field t-2 g -e g 
splitting, orbital differentiation happens once entered the 
Janus regime, where the orbitals closest to half-filling are 
more correlated (2^ . 
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